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Abstract: We study bases of Pfaffian systems for A-hypergeometric system. 
Grobner deformations give bases. For hypergeometric system associated to a 
class of order polytopes, these bases have a combinatorial description. These 
theoretical study gives us algorithms of transforming a given ^-hypergeometric 
system into a Pfaffian system. The size of the Pfaffian system associated to a 
subclass of the order polytopes has the growth rate of the polynomial order. 

1 Introduction 

Let be a vector valued function in xi , . . . , x„ . We suppose the length of F is 
r. Let Pi{x), i — . . . ,n he r X r matrices satisfying 

dxj ^ dxi ^ 

for all i ^ j. The system of linear differential equations 

dF 

— = Pi{x)F, i = 1, . . . ,n 

OXi 

is called a Pfaffian system. We also call the system of the linear differetial 
operators ^ — Pi a, Pfaffian system. The number r is called the size of the 
Pfaffian system. For a given zero dimensional left ideal in the ring of differential 
operators with rational coefficients, it is well-known that an associated Pfaffian 
system can be obtained by a Grobner basis method (see, e.g., [121 Appendix]). 
A geometric method to find a Pfaffian system associated to a given definite 
integral with parameters is the use of twisted cohomology groups (see, e.g., [2])- 
In |12| , new methods in statistics are proposed. They are holonomic gradient 
method (HGM) and holonomic gradient descent (HGD). The HGM is a method 
to evaluate numerically the normalization constant, which is a function of pa- 
rameters X, for a given unnormalized probability distribution, and the HGD is a 
method to make a maximal likelihood estimate by utilizing the HGM. The key 
step of these methods is to construct a Pfaffian system associated to the nor- 
malization constant. The size of the Pfaffian system determines the complexity 
of the HGM and the HGD (see, e.g., [TO]). 
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^-hypergeometric systems have been studied for the past 20 years (see, e.g., 
[5], [6], [T5|) with several apphcations in various fields. We will add a new 
motivation to study it; we note that it can be regarded as a system of differ- 
ential equations for normalization constants of a class of exponential family of 
probability distributions (see, e.g., ((TU)) ). Following these studies and the new 
motivation, we pose the following questions to apply A-hypergeometric systems 
to statistics. 

1. Give an efficient method to find a Pfaffian system associated to a given 
A-hypcrgcomctric system. 

2. Find a nice subclass of ^-hypergeometric systems of which associated 
Pfaffian systems have a moderate size. 

We will give an answer to these questions in this paper. Let us be more precise. 
Let Rn be the ring of differential operators with rational function coefficients 
in n-variablcs. The first step to find a Pfaffian system associated to a given 
left ideal / in i?„ is to give a basis of i?„// as a C(x)-vector space. When the 
set {ui, . . . , Ur} is the basis, there exists a matrix Pi{x) such that d/dxiU = 
PiU mod/, U ~ (ui, . . . , UrY' . We can show that d/dxi — Pi is a Pfaffian system. 
We call . . . , Ur} the basis of the Pfaffian system. We study bases of Pfaffian 
systems for A-hypergeometric systems. We will show that Grobner deformations 
give bases (Theorem [T|) . For hypergeometric system associated to a class of 
order polytopes (see, e.g., [9]), these bases have a combinatorial description 
(Theorem [5]). These theoretical study gives us algorithms of transforming a 
given A-hypergeometric system into a Pfaffian system (Section |6]). 

The size of the Pfaffian system associated to a subclass of the order polytopes 
has the growth rate of the polynomial order (Theorem 15]). We expect that 
this together with our algorithm yields a new class of an exponential family of 
probability distributions for which we can apply the holonomic gradient method 
(HGM) and the holonomic gradient descent (HGD) efficiently. 

2 Bases for the Pfaffian System 

Wc denote by A = (oy) &dx n-matrix whose elements are integers. We suppose 
that the set of the column vectors of A spans Z''. Let si, . . . , be indetermi- 
nates. Following notations in [TS], we denote by i/A[s] a left ideal generated in 
the Weyl algebra 

D[s\ = C[si, . . . ,Sd]{xi, . . . ,Xn,di, . . . ,dn), di = d/dxi 
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by 



n 



7 , '^ij^j^j 



- Si 



(1) 



n 



71 



(2) 



(with u,v G Ng running over all u,v such that Au 



Av). 



Here, Nq = {0, 1,2,.. .}. We call the ideal generated in C[di, . . . , 9„] by the 
elements of the form ([2]) the affine toric ideal and denote it by I a ■ We denote 
by Ei — Si the operator ([T]). For complex parameters /3i, the system of linear 
differential equations {Ei ~ Pi)f = (i = (9" - 9")/ = {Au = 

Av) is called the ^-hypcrgcomctric system of differential equations or shortly 
the A-hypergcomctric system. We also sometimes call the ideal the A- 

hypergeometric system (with indefinite parameters). 

Let Rn be the ring of differential operators with rational function coefficients 



We are interested in bases of Rn/{RnHA[s\) as the vector space over the field 
C{s,x). Any basis of the vector space yields an associated Pfaffian system 
or an integrablc connection associated to Ha [s\ ■ Let ui , . . . , be a basis of 
Rn / {RnH a\s]) ■ For Uj, there exist rational functions Pi^ S C{s,x) such that 
diUj = Y^\^-^^p^jUkU\od RnH a[s\- The system of differential equations diF = 
{Pij I 1 < Ji ^ f )F where F is a vector valued function of size r is called a 
Pfaffian system and {ui} is called a basis of the Pfaffian system. 

Bases can be described by those of simpler quotients, of which denomina- 
tor ideals are nothing but Grobner deformations of -ffyi[s], as in the following 
theorem. 

Theorem 1 Let w e Z" be a generic weight vector for the affine toric ideal I a 
such that Aeghitu{I a) =deg/^. Letui, . . . ,Ur be a monomial basis of Rn / {Rn J) 
where the left ideal J is generated by u1tu{Ia) and Ei — Si, i — 1, . . . ^ d in Rn- 
Then, the set {ui, . . . , u, } is a basis of the vector space Rn/{RnHA[s])- 

Proof. We denote by r the normalized volume of A. Since Si are inde- 
terminates, the holonomic rank of J and i?yi[s] are r by Adolphson's theorem 
(see, e.g., [T], |15|). In other words, we have dimc(^s .^.^Rn/ {RnH a[s]) = r and 
<iimc(s,x)Rn/J = r. 

We may assume that Ui are expressed as a monomial in terms of Euler 
operators 6j ~ Xjdj. When we regard J as a system of linear differential 
equations, it has r linearly independent solutions of the form where p € 
C(s)". We denote them by gi = a;'''^*^ i = l,...,r. Since gi are linearly 
independent solutions, the Wronskian determinant det(Mi ^gj) is not identically 
equal to 0. The solution gj can be extended to a solution fj of such 



C{s,x){di,. . . , 



dn). 



(3) 
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that gj is the leading monomial of fj with respect to the weight vector w (see, 
e.g., [15l Chapters 2 and 3]). The series fj is expressed as fj = gj J^teM- C'fa;^, 
Co = 1 where Mj denotes the set of the lattice points in a cone and Ce is 
a constant belonging to C(s). The series converge in the space of convergent 
power series gj ■ 0{U){Mj} where U is an open set in the s-space and 0{U) is 
the space of holomorphic functions on U [13]. We replace Xi by Xif^' for all i 
in fj and denote by xt^ the vector {xif^^, . . . jXnf^")- From the construction 
algorithm of fj, we may assume that fj{xV") = gj{xt^){l + 0{t)) when < — >• 
as a function of t when x is fixed and s lies in U . 

Let us prove W = det(Mi • fj) ^ 0. We denote by Ui{p{j)) the constant 
[ui • x^^^^) / x''^^\ Under this notation, we have x'^^^^Ui • gj = Ui{p{j)) and 

{x-'>^'\u, . = u,{p{j)+£)C,xV^ (4) 

Note that > Q for (. ^ Q and £ S Mj . Therefore, we have 

detix-P'^^^u, • fj){xt'") = dct(a;-''(^')u, • gj){xr') + 0{t) (5) 

from (|4|) when a; is fixed and t — > 0. This implies that the Wronskian determinant 
Aet{ui»fj) = (^rij det(a;^''(^'ui»/,) is not identically equal to 0. Therefore 

Ui arc linearly independent in Rn / {RnH a[s\) ■ Q.E.D. 

Remark 1 We may expect that when the set {ui\ is a set of the standard 
monomials of J, then it is also a set of the standard monomials for i?„7?A[s]- 
However, it is not true in general. In fact, if the set {ui} is the set of the 
standard monomials of J for a term order >- in Rn, then there exist rational 
functions such that 9" — J2j '^a'^j '= RuHaIs] for any 9" ^ {ui}, which can 
be a Grobner basis of RnHA[s] for a weight vector in the Grobner cone defined 
by a ■ w > w ■ (exponent of Uj) for a^^ ^ if the cone is not empty. 

Let M be a monomial ideal in C[d]. When M is generated by d", the 
distraction M C C[9] is generated by n"=i — 1) • ■ • (0i — + 1) where 
9i = Xidi [ini p. 68]. Put M = v[\w{Ia). Then, the ideal J in the Theorem[T]is 
generated by M and ^ij^j ^ * = 1, . . . ,d [HI Sec. 2.3, Prop. 3.1.5]. 

This leads us to the following corollary. 

Corollary 1 Retain the assumptions of the Theorem [H The set of the mono- 
mial basis of C{s)[6]/ J , where J is the ideal generated by M and ^ciijOj — Si, 
i = I, . . . ,d in the polynomial ring C{s)[9], gives a basis of Rn/ RuHaIs] by the 
replacement 9i = Xidi. 

Wc denote by (0,1) the weight vector such that Xi has the weight and di 
has the weight 1. The initial form in(o,i)(^) is denoted by a{l) following the 
standard notation of the principal symbol in the analysis. a{di) is set to ^i. For 
D[s) = C{s){x,d), we have a{D{s)) = C{s)[x,£_]. The induced map from i?„ 
to C(s,a;)[^] is also denoted by a. 
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Theorem 2 ([5]) Suppose that I a is homogeneous and Cohen- Macaulay (,i.e., 
C[d]/lA is Cohen- Macaulay). Then, the ideal a{D{s)HA[s\) is generated by 

n 

(j{Ia) and a^Ei — Si) = '^ij^j^.j: i — 1, ■ ■ ■ , d. 

We assume that Ia is Cohen-Macaulay in the remaining part of this section. 
Suppose that the set of monomials pi, . . . in ^ is a basis of C(s, x)[S,]/a{D{s)HA[s]) 
as a vector space over C{s,x). The following proposition gives a different way 
to give a basis of R„/RnHA[s]. 

Proposition 1 The set of operators pi = p^^^Q, i = l,...,r is a basis of 
Rn/ RuHaIs] as a vector space over C{s,x). 

Proof. We note that the dimension of Rn/ RuHaIs] as the vector space and 
the dimension of C{s,x)[^]/(j{D{s)Ha[s\) agree and are r. Suppose that it is 
not a basis. Then, there exist non-zero polynomials Ci{s,x) 6 C(s,x) such that 
7^ X^CiPi e RuHaIs]. By using the relation i?„i/^[s] n D{s) = D{s)Ha[s\ 
(Theorem |4]), wc may assume that q S C[s,x] and '^ciPi G iJ^is]- Take the 
principal symbol of it, then we have 

^C^C^p^) = ^ C^Pi{£,)-.p. 

i such that (0, l)-order of pi is maximal 

By the Theorem [51 we have p S a{D{s)HA[s]). It contradicts to that pi, i = 
1, . . . , r is a basis. Q.E.D. 

Wc denote the toric ideal (^" — | Au = Av) in the ^-spacc also by I a as long 
as no confusion arises. A natural question coming from the discussions above is 
a relation between bases of C(s, x)[^]/C(s, x){Ia, ^ij^j^j | * = 1, • ■ • , and 

of the simpler quotient C(s)[6']/C(s)[6'](inu,(/^), ^ o-ij^^j — | i = 1, . . . , d). 

Giving the weight 1 to diS, we define Fk — {}2,\a\<k^a{x)d°' | e K{x)} 
where K is C(s). Let / be a left ideal in i?„. The map a gives an isomorphism 
between i^-vector spaces 

a : Fk/il n Fk + F^^i) n a{Fk) + aiFk-i)) 

for fc = 0, 1, 2, . . .. Note that the set of the standard monomials of Rn/I with 
respect to an order >-i gives a basis of gr(i?„//) = ®^^qF}~/{I n + Fk-i). 
Let 7 be a map which translate a basis {ui\ in the Theorem [1] to a basis of 
^{Rn/ RuHaIs]). Note that when UiS are the standard monomials of RuHaIs] 
with respect to an order we have 7 — id. Then, we have a C(s, a;)-linear 
isomorphism 

a o 7 : Rr,l[RnHA[s\) ^ gr(i?„/i?„i7^[s]) ^ gv{G{Rn) / a{RnHA[s\)). (6) 

By using the map 7, we obtain the following theorem in the commutative alge- 
bra. 
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Theorem 3 Suppose that the affine toric ideal is homogeneous and Cohen- 
Macaulay. Let w he a generic weight vector for the toric ideal I a ■ Put M = 
(Ia)- 

1. Let M he the distraction of M which lies in the commutative ring C{s)[6]. 
Then, a hasis of the quotient space of C{s)[9] hy the ideal generated hy M 
and cLijOj — Si, i ~ 1, . . . ,d gives a hasis of C(a;)[^]/ {I a, Y1 '^ij^j^j I * = 
1, . . . , d) as the C{s) -vector space hy taking the image hy a o ^ . 

2. Lf M = 'niw{lA) is square free, then a hasis of the quotient space o/C[^] hy 
the ideal generated hy M and the parameter system ^ o-ij^j, i = 1, ■ . ■ ,d 
gives a hasis of C{x)[^]/ {La, '^ij^j^j \ i = ^, ■ ■ ■ , d) hy taking the image 
hy a o J. 

Proof. (1) Applying the Corollary [1] and the Theorem [2l wc obtain the 
conclusion by the definition (|6]). (2) When M is square free, M agrees with M. 
Then, (2) follows from (1). Q.E.D. 



3 Ha[s] as the Annihilating Ideal of Hypergeo- 
metric Series 

We denote by {—w,w,0) the weight vector in D[s] such that the weight —Wi 
stands for Xi, Wi for di, and for all s^'s. Let -< be a block lexicographic order 
such that si y S2 y ■ ■ ■ y Sd >- other variables. The order for "other variables" 
may be any term order. 

Proposition 2 (\15[ Theorem 3.1.3]) Let Ga he a Grohner hasis for the homo- 
geneous toric ideal La for the order -<(-w.w.o)- Then the set 

GaLiIs,^ E,\i^l,...,d} 

is a Grohner basis of L[a[s] in D[s] for the order '<(~w,w,o)- 

A proof of this proposition can be done by the Buchberger algorithm in the 
homogenized Weyl algebra D[h,s] where the weight for the homogcnization of 
Si is 2 and the weight for other variables is 1. 

Let w be a generic weight. When La is homogeneous, there exist r = vol(A) 
independent solutions of the form x''* i = l,...,r for the initial system 
in(-u,,u,,o)(^A[s]) = (in„(/yi),Si - Ei\i ^ I, . . . ,d). 

Proposition 3 Let d £ D[s] he a differential operator which annihilates all 
xP'^'l Then, e e in(_^,^,o)(^A[s]). 

Proof. Using Ej — Sj, we can delete the variable Sj from £. Then, we may 
assume that the operator £ does not contain the variable s. We decompose 
the operator £ into a sum of Euler operators as £ = where 9j = 
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Xjdj. Since we have I • x^'^'^'^ = ^ ia[p{s))x°'~^''^'^^ = 0, we may assume that 
la{0) • xP'^'^'> = for any i and a. This leads a system of algebraic equations 
^a{pi{s)) = OinC[s]. We can show that £a (0) & D ■m.^^A)- Since -D-in„(/A)n 
C[0](the distraction) + [AO — s) is a radical ideal of which number of solutions 
is r [T51 Lemma 3.2.4, Theorem 3.2.10], we obtain the conclusion. Q.E.D. 

The solution of the initial system can be extended to a solution of 

Ha[s\. We denote the solution by fi. The following theorem is shown by 
L. Matusevich [TT]. We give a slightly different proof. 

Theorem 4 (Ul'l) Let L he a differential operator in D[s\. If L • fi — Q for 

any i G {1, 2, . . . , r}, then L G Ha[s]. 

Proof. We express the operator L as L = ^i^s". La G D. Using Ej — Sj, 
we can delete Sj from L. Therefore, we may assume L € D. The operator L 
can be expressed as ^x^S*. Let Ga be a Grobner basis of I a for the generic 
weight vector w. By reducing each 9'? by the Grobner basis Ga, we may assume 
that 

L= J2 

9'? is a standard monomial 

Suppose that L ^ 0. Since i • = 0, we have 

in(-t„,t„,o)(i) • Startii,(/,;) = 0. 

We note that Starttu(/,;) — xP^^'^\ We apply the Proposition [3l Then, we 
have 'va-(-w,w.o){L) G '™.(-w,wfl){HA[s\)- Let >- be a term order such that si >- 
■ ■ ■ y Sd >~ (other variables). The monomial operators ui.uj{Ga) and the Euler 
operators Sj — Ej are a Grobner basis for the order y(-w,w,o) (Prop. [2]). Since, 
i'n{~w,w,o)iL) is irreducible by this Grobner basis, we conclude that L ~ 0. 
Q.E.D. 



4 74-hypergeometric Systems for Order Polytopes 

First, recall what the order polytope of a finite partially ordered set (|9l p. 115]) 
is. Let P ~ {ai, . . . , ad} be a finite partially ordered set with \P\ ~ n. A poset 
ideal of P is a subset a of P such that if a G a and h £ P together with b < a, 
then b € a. Thus in particular the empty set as well as P itself is a poset ideal of 
P. Let J{P) denote the distributive lattice (y, p. 118]) consisting of all poset 
ideals of P, ordered by inclusion. For Example, if P is the disjoint union of 
two chains of length 2 and length 3 drawn in Figure [I] then L ~ J{P) is the 
distributive lattice drawn in Figure [51 

Let ei, . . . ,ed denote the standard unit coordinate vectors of R''. If /3 is a 
subset of P, then we write Wfj for the (0, l)-vector g/9 ^« ^ ^^'^ order 
polytope 0[P) C R'' of P is the convex hull of the finite set {wa '■ a G J{P)}. 
Its dimension is dimC'(P) = d. 
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Figure 1: 



Figure 2: 



Let K = C{{£,a}a£j{P)) denote the rational function field in |J'(P)| vari- 
ables over C. Let A = K[ti, . . . ,td, s] denote the polynomial ring in n variables 
over K. If /3 is a subset of P, then we write up for the squarefree monomial 
OaiG/s ^[^(P)] denote the subalgebra of A which is generated by those 

squarefree monomials up with (3 £ J{P)- The semigroup ring K[0{P)\ is in- 
troduced by [8]. We caU K[0{P)] the toric ring of 0{P). The KruU-dimension 
oiK[0{P)] is d+ 1. 

Let K[{xa}a^j[p)] denote the polynomial ring in \J{P)\ variables over K 
and define the surjective ring homomorphism tt : K[{xa}aeJ[P)\ ^ ^^['•^{P)] 
by setting 'K{xa) = Ua- Its kernel Io{P) is called the toric ideal of 0{P). It is 
known [8] that Io{P) is generated by those quadratic binomials 

XaXp — Xa/\pXa\/p (7) 

such that a and /3 are incomparable in the distributive lattice J{P)- We fix 
an ordering < of variables of K[{xa} a£j(P)\ with the property that if a > 
P in J^{P), then Xa < xp. Let <rcv be the reverse lexicographic order on 
K[{xa}a£j(P)] induced by the ordering <. In [8] it is shown that the set of 
binomials ([7|) is the reduced Grobner basis of Io{p) with respect to <rev Thus 
ii^<,ov (^o(P)) is generated by those squarefree quadratic monomial XaXp such 
that a and /3 are incomparable in J{P). 
Let 

6, = ^ (,aXa, 1 <i <d 

together with 

aeJiP) 

It follows that the sequence (6'o, ^i, • ■ • , dd) is a system of parameters of both 
the residue rings K[{xa}a(^j{p)]/ Io{P) and K[{xa}a£JiP)]/'^n<,^^{Io(P))- A 
fundamental question is to find the standard monomials of the 0-dimensional 
residue ring 

^[{a^a}aej(p)]/(in<_(/o(P)), 00, ^i, • ■ • , Oa). (8) 
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In general, however, the question is very hard. When P is the disjoint union of 
two chains, a complete answer can be given. 

Let Pp q denote the disjoint union of two chains Cp : ai < • ■ • < Op of length 
p ~ 1 and Cq : bi < ■ ■ ■ < bq oi length q — 1. Let a^.j, where < i < p and 
J 5; be the posct ideal {ai, . . . , a^, &i, . . . , 6^}. In particular ao.o = 0- 
Write for , and x^.j for Xa^ j ■ Let 

Ot* = ^ ^k.jXkj, <i <P 

i<k<p, a<j<q 

and 

0<i<p, j<l<q 

Thus in particular 

0<i<p, 0<j<q 

Let A"[x] = K[{xi,j}a<i<pfi<j<q]. Then the residue ring (|5]) of Pp^g is A'[x]/Jp 
where 

} i<k, e<j, {^i*}o<i<p, {(^*j}o<j<q)- 

Let <rev denote the reverse lexicographic order on K[x\ induced by the ordering 
of the variables 

a;0,0 > 2:1,0 > 2:04 > 2:2,0 > 2:1,1 > Xo,2 > • • • > Xp^q. 



Lemma 1 In /i'[x]/in<_.^^ ( Jp.g) ~ together with 

Xz,o = xoj = 0. 

Proof. Let i < i'. Then 

0i'*Xij-Xij{ ^ ^kjxk,e) 

i'<k.j<e 

belongs to {{xi^.jXk,i}i<k,e<3)- Hence 

i'<k,j<l 

belongs to J. Thus its initial monomial XijXi'j belongs to in<j,,^ ( Jp.g). Let 
i = i'. The initial monomial of 

Gi^Xi^i ^i^i ^^i ^ ^ ^i.kXi^k 
k=0 

is x^j. Hence j £ in<^.^^ ( Jp,,). Similarly, XijXij' S in<^^^ ( Jp,,). 

Moreover, since and belong to J, their initial monomials Xi,o and ccg j 
belong to in<_.^^(Jp, g). Q.E.D. 
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Let Sp^q denote the set of squarcfrec monomials of K[x.] of the form 
where 

< ii < i2 < ■ ■ ■ < ir < P, < ji < j2 < ■ ■ ■ < ir < q, r = 0, 1, 2, . . . . 

Theorem 5 The set of standard monomials o/ in<^^^ ( Jp q) is equal to Sp_q. 

Proof. Lemma [T] guarantees that each standard monomial must belong to 
Sp^q. In [3] it is proved that the number of standard monomials of degree r 
coincides with the number of maximal chains of J^{P) with r descents. Recall 
that the descents of a maximal chain 

of J{P) are those ai^^j^ with 1 < fc < p + g such that 

T-k-i = ik < ik+i, jk-i < jk = jk+i- 

Now, given a squarefree monomial ([9|) of degree r, we can associate a unique 
maximal chain whose descents are 

in the obvious way. (See Figure [S]) Hence the number of squarefree monomials 
^ of degree r is less than or equal to that of standard monomials of degree 
r. It then follows that Sp_q is the set of standard monomials of in<^^^ ( Jp,g), as 
desired. Q.E.D. 



Remark 2 Let so,Si*,s*j,, where I < i < p and 1 < j < q, be indeterminates 
over K = C{{S,a}aeJ(P)) K' = A'(so, {s«}i<i<p, {s,j}i<j<g). Let = 
0i* — Si-f and = 0fj — Sitj, where 1 < i < p and I < j < q, together with 
^0* = ^0* ^ •^0 and = — sq . Then the proof of Lemma Q] as well as that 
of Theorem [5| is valid for the ideal J'^ ^ of K' [x] , where 

J'p.q = {{Xi,jXk,l>}i<k,t<3-,{Q'i^}vi<i<p-,{Q't.j]o<j<q)-, 

and the set of standard monomials of m^^^^{Jp q) is also equal to Sp_q. 

Remark 3 For Jp^q, p = 1, the map 7 in (d^ can be the identity map. Then, 
Si^q gives a K-hasis of the vector space 

K[{xa}aeJ(P)\/{Io(P),So, 6*1, • • • , Od) 

where P = Pi^q. Let us prove it. We have Si^q = {l,Xi^i,Xi^2, ■ ■ ■ ,xi_q} then, 
by the CorollaryUl {1, di,i, 9i.2, . . . , ^i.^} is a basis of the Pfaffian system. For 
any d" , of which order is larger than 1, there exist a" G C(s, a;) such that d°' — 
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Figure 3: 



Oq — X]j=i '^f'^ij 6 RuHaIs]. There are also operators of the forms xqjOqj + 
xijdij — (a constant), a;o,iC^o,i + X]J=i ^i.j^ij ~ constant) in i?„i7^[s]. T/iey 
are a Grobner basis in Rn with the (graded) reverse lexicographic order The 
leading terms are underlined. By the Remark\^and the Theorem\^ we obtain 
the conclusion. 

Wo now turn to the discussion of the normalized volume of order polytopcs. 
It follows from [17] that the normalized volume of the order polytope 0{P) is 
equal to e(P), the number of linear extensions of P. Recall that an antichain 
of P is a subset B oi P such that if a and b belong to B with a ^ b, then a 
and b are incomparable in P. The width of P is the supremum of cardinalities 
of antichains of P. The length of a chain C is |C| — 1. The rank of P is the 
supremum of lengths of chains of P. 

Lemma 2 Fix positive integers d and r. Let P be the disjoint union of d chains 
Ci, . . . ,Cd and suppose that the length of each chain Ci with 1 < i < d is at 
most ?■ — 1. Then there exists a polynomial f{n) in n of degree r{d — 1) such 
that e{P) is at most f{n), where n = \P\. 

Proof. Let £i denote the length of C;. Then the number of linear extensions 
of P is 

/ n \ n! 
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Since id = n — X^iLi Ci ^ n — r{d — 1), it follows that 



' ' - £d! - (n-r(d-l))!- 

Let 

f{n) ^ n{n - l){n ~ 2) ■ ■ ■ {n - r{d - 1) + 1), 

which is a polynomial in n of degree r{d — 1). Then e(P) < /(f^), as required. 
Q.E.D. 

Theorem 6 Fix positive integers d and r. Let P be a finite partially ordered 
set and suppose that there exists a chain C of P such that 

(i) the width of P \ C is at most d — I; 

(ii) the rank of P\C is at most r — 1 . 

Then there exists a polynomial f{n) in n of degree r{d ~ 1) such that e(P) is at 
most f{n), where n ~ \P\. 

Proof . Since the width of P \ C is at most d — 1, Dilworth's theorem [4] 
guarantees the existence of d— 1 chains Ci, . . . , Cd-i of P\C, where the length of 
each d is at most r - 1, such that P\C = C1UC2U ■ ■ -U Cd-i and Ct n Q = 
for i 7^ j. Hence there exists a partially ordered set Q which is the disjoint 
union of d chains C(, . . . , C^, where the length of each with 1 < i < is at 
most r — 1, such that there is an order preserving bijection (p : Q P. Hence 
e{P) < e{Q). Thus the desired result follows from Lemma[5J Q.E.D. 

Let us come back to our discussion on the A-hypergeometric system. We note 
that we have used a standard convention of variable names in combinatorics and 
commutative algebra, which is not compatible with that in A-hypergeometric 
systems and D-modules. Here is a table of the correspondence. 



Combinatorics 


^-hypergeometric system 


X 


^ or 9 


e 


X 


6 stands for a parameter system. 


9 stands for the Euler operator 



The results of this section together with those in the Section[2]yield the following 
claims on Pfaffian systems of ^-hypergeometric systems. 

1. By Theorem [5] and Corollary [1] the set Sp^q\xij-^dij gives a basis of the 
Pfaffian system for the A-hypergeometric system defined by the order 
polytopc 0{Pp^q), which agrees with the Aomoto-Gel'fand system E{p + 
l,p+l + q+l) i2J. 

2. By Theorem [6l we can regard the rank e{P) of the hypergeometric sys- 
tem associated to the order polytope 0{P) (n = \P\) has the polynomial 
growth property with respect to n. 
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5 Bases of Twisted Cohomology Groups 



Let Ai = (ai, . . . ,a„J, . ..,Ah = (a„j^_i+i, . . . ,a„J, a,, G Z'". To each matrix 
Aj we associate a generic polynomial fj{x,t) = X]i=nj_i+i ■ 1^°^ parame- 
ters ai, . . . , afe and 71, ... , 7^, we consider the integral 

$(a,7;x)= / P{x,t)dti---dt^, P{x,t) fj{x,tTH^ (10) 

for a suitable twisted cycle C. The function $ is satisfied by the A-hypergeometric 
system for 



/ 1 





A 



\ ai 



1 



\ 






and /3 = (ai, . . . , a^, —71 — 1, ... , —7m — 1)"^ [6]. We regard P as a function in 
t — {ti, . . . , tm) with the parameter vector x. Put D{P) = P\a=^=i- Define the 
connection V with rational function coefficients by 



V 



dti 



(11) 



where d is the exterior derivative with respect to the variables ti, . . . ,tm- When 
fj=Yl Xif^^^s are linear polynomials in a generic position, it is well-known that 
the space of the cohomology group H™{n* {*D{P)),V) is isomorphic to the 
germ of the solution space of the ^-hypergeometric system when the parameter 
x lies in a domain out of the singular locus and 01^,7^ £ C \ Z. This isomorphism 
gives the following Theorem from our Theorems [T] and [5l 

Theorem 7 Let A be the matrix standing for the order polytope 0{Pp^q) and 
he the set of standard monomials given in Theorem\^ Then, the 
set of the rational forms 



Ui • P 



P 



-dti ■ ■ ■ dtm, i = 1, . . . , r 



(12) 



is a basis of the twisted cohomology group i7™(fi*(*D(P)), V) when ai,^j G 
C\Z. 

This theorem is a different presentation of the celebrated work of K. Aomoto who 
gave a basis of the twisted cohomology group associated to a generic hyperplane 
arrangement (see, e.g., [H Theorem 9.6.2]). 2012/ 
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6 Algorithms 



Let / be a zero dimensional left ideal in the ring of differential operators with ra- 
tional function coefficients, which we denote by i?„ = Kiiji, . . . , yn){di, . . . , 9„), 
or Rn{K), where K is a. field of characteristic and di = d/dyt. We also denote 
the ring of differential operators in the cc-variable by i?„ as in the previous sec- 
tions as long as no confusion arises. We assume that the holonomic rank of / is 
r. We denote by = (ui, . . . , UrY' ^ column vector of which elements form a 
basis of the vector space i?„/ J over the field K{y) = K{yi, . . . , y^) {F is not a 
vector valued function in this section). A Pfaffian system of factored form [10] 
for / and F is Qidi — Pi, i = 1, . . . , n satisfying 

QAF = P,F vaodJ (13) 

where Qi ^md Pi are r x r matrix with polynomial entries and Qi is invertible. 
The operator di in the definition may be replaced by Oi = yidi. Let G be a 
Grobner basis of / and F a vector of which elements form a set of the standard 
monomials of G. Then, the matrix Q~^Pi can be obtained by computing the 
normal form of diF by G. This matrix is called the companion matrix in the 
theory of Grobner basis. However, the computation of Grobner basis in i?„ is 
very heavy task and we can derive a Pfaffian system only for relatively small 
sized /. 

Let us consider the problem of transforming an A-hypcrgcomctric system 
into a Pfaffian system. The translation is possible by computing Grobner basis 
of RnHA[s] as we have explained. We propose two efficient methods which work 
for A-hypergeometric systems. 

The first method is to utilize the rank theorem and the Hilbert driven algo- 
rithm [18j in the simplest form. 

Algorithm 1 

Input: a d X n matrix A with integer entries such that TiA 
1, . . . , d) is a free 'L-module of rank d. 
Output: a Pfaffian system. 

1. Evaluate the rank r of Ha[s\ by the normalized volume of A. 

r = the normalized volume of A. 

2. Apply the Buchberger algorithm for Ha [s\ in i?„ . We stop the computation 
when the Hilbert function is equal to r. Let the output be G. 

3. Construct the Pfaffian system di — Q^^Pi by the normal form algorithm 
in Rn by G. 

The set G is a Grobner basis since the holonomic rank of RuHaIs] is r by the 
Adolphson's theorem [1]. Therefore this algorithm outputs the correct answer. 



Ej=i Z(a»j I i = 
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Figure 4: Ci 



Example 1 When the point configuration is defined by the lattice in the Figure 
m the matrix A (the vertices of the order polytope associated to the figure) is 

1\ 

1 

1 ■ 
1 / 

Let us derive a Grobner basis of R„Ha[s] and Pfaffian systems by the Algorithm 
[TJ The Hilbert driven Buchberger algorithm takes 8.22s and it takes 168.86s to 
obtain a Pfaffian system (for all directions di, . . . ,ds) from the Grobner basis. 
If we do not use the Hilbert driven method, the Buchberger algorithm takes 
3457.56s. It is 421 times slower than our algorithm. (The timing data are 
taken on a machine with the AMD64 CPU (3.33GIIz) and 32GB memory, the 
computer algebra system Risa/Asir and the package yang.rr [14].) 

The second method is an algorithm of deriving a Pfaffian system without 
computing a Grobner basis of the ^-hypergeometric ideal J = RnHA[s] where 
K = C(s). We use series solutions of the ^-hypergeometric system to construct 
a Pfaffian system. The correctness of our method follows from the discussions 
in the sections [2] and [3l We list facts on which are used in our algorithm. 

1. A {—w, w, 0)-Gr6bncr basis of the hypergeonietric ideal i?yi[s] in D[s] con- 
sists of the Euler operators Ei — Si and a w-Grobner basis of the toric ideal 
Ia [m Theorem 3.1.3]. Here, —w is the weight for x and w for 9, and 
for the variable s (Prop. [2]). 

2. If we have a Grobner basis of Ia, then we can efficiently construct se- 
ries solutions of i?yi[s] (see, e.g., [15j Theorem 2.5.14, Section 3.4] or [71 
Theorem 6.12.14]). 

Let / be a series solution of i?yi[s]. We may assume that / belongs in 
xp'K[[C* n Z"]] where K = C{s), p' e K and C is a Grobner cone defined by a 
Grobner basis of Ia [HI Theorem 2.5.14]. We construct a unimodular cone U 



Xi X2 X3 X4 X5 Xq Xj 



A = 
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1 


1 


1 


1 
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1 





1 
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1 











1 





1 





1 











1 





1 
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contained in the Grobner cone C by the algorithmic procedure for the resolution 
of singularities of toric varieties (see, e.g., [16]). The dual U* is also a unimodular 
cone. Let m{i) be rays of U* such that U* n Z" = ^Z>om(«). The matrix 
(m(l), . . . , m(ri,)) is a unimodular matrix. By a unimodular transformation 
of monomials yi — x"^^^\ we may assume that / in terms of y belongs to 
yPK[[yi,...,y„]]. 

Example 2 Using the resolution algorithm should be the final strategy after 
all other methods fail. Let us find a unimodular cone in a Grobner cone for the 
A in the Example [T] without the resolution algorithm. The Grobner basis for 
the reverse lexicographic order is 



{X2X3 



2:4X5 

and the Grobner cone is 



XiXe,X2X7 
XiX3,X5Xe 



■ XiXs,X3X4 
X2Xs,X5X7 



XiX7,X3Xg 

X3Xs,xeXY 



XlXs, 
XiXg) 



c 



It) G 



Cone 



Wc define Ci {i 
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1, . . . , 9) by each vector appearing above expression from left 
to right. In this case, we do not need to apply the resolution algorithm and 
have a simpler method. In fact, the cone (cq, ci, C2, C3, C4, C5) gives an affine 

(1,0)^, 



toric ideal (aoaia2 ^a3a4ac,) where cq 



and a," arc variables 



corresponding to q. It has the unimodular decomposition because it has a 
square-free initial term. Matrices (ci)i<i<g \ cj {j = 1, 2, 3, 4, 5) are unimodular 
and the cone C contains cones generated by them. In the case of j = 1 , wc can 
take 



U = Conc((ci)i<i<9 \ ci) Cone 



/ -1 
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and its dual cone ?7* is 



U* = Cone 



V 













-1 



1 1 



-1 




1 
1 



-1 / 



Let Y'^ be an ideal in C[?;] generated by the monomials of degree k. In the 
sequel, we assume that a Pfaffian system is expressed in terms of the variable 
y. We assume the elements of F are monomials of Euler operators. We denote 
by F{f) a vector obtained by acting operators in F to the series /. 6*^ is the 
operator in the y variable standing for Xid/dxi. Then, we have 



(14) 



Since / G y''C [[?/]], any element of O'^Fif) and F{f) belongs to yPK[[y]]. There- 
fore, wc have 



Q^iy-^e'^Fif) modY") = P^iy-PF{f) modY'^) 



(15) 



in the polynomial ring K[y] modulo Y'^ for any natural number k. Moreover, 
multiplying a suitable polynomial c(s) in s to F{f), we may assume that it 
holds in C[s,y]. In summary. Pi and Qi appear as a syzygy among truncated 
solutions as in ([T5|). This leads us the following algorithm to obtain Qi and Pi. 



Algorithm 2 



1. Fix a number k^. Determine a vector F of monomials of Euler operators, 
which is a basis, by the Corollary [1\ (or other methods in sections [H and 

2. Solve the polynomial syzygy equations 

q,-c{s){y~Pe[F{f) modF'^) - p, ■ c{s){y-P F{f) modr'=) 
+ ^'^.^ ■ = (16) 

\a\=k 

in C[s,y] for unknowns qi, pi, r^.a where qi,Pi are polynomial row vectors 
of which length are r = rank(iJ^[s]) and f runs in the space of solutions 
of Ha[s] and k runs from 1 to kij. Here, F(f) is expressed in the variable 
y. Note that c(s) and p depends on f and k and are chosen so that 
c{s){y-Pe[F{f) modr'=) and c{s){y-P F{f) modF*^) belong to C[s,y]. 
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3. Construct matrices Qi and Pi from the syzygies. If we cannot construct 
them or Qi is not invertihle, then increase fco and go to 2. 

4- Apply the candidate Qi9[ — Pi to the series solution F{f) To.odY^'^^ trun- 
cated at the degree k ^ N + \ > fco • Determine remaining undetermined 
constants by setting {Qi9[ — Pi) • {F{f) niodF'^) = mody'"'. // there is 
no solution, then increase fcg and go to 2. 

5. (Check the candidate.) Reduce Qi9[F — PiF expressed in the variable x by 
a {—w,w,0)-Grobner basis of Ha[s] by the normal form algorithm. The 
details of this step will be discussed in the algorithmic If the result is yes, 
then we are done else, increase fcg ciiT-d- go to 2. 

Remark 4 The operator 6[ in the Algorithm may be replaced by any operator 
of the form ^ y'^pa{9y) where Pa is a polynomial in 9y-^ , ■ • • , 0y„ to find a "con- 
nection" to the direction "J^V^Pai^y)"- This means that the Algorithmic can 
also be used to express the gradient grad ft and the Hessian Hess ft in terms of 



We note that we may solve the syzygy equation by the Grobner basis in the 
ring of polynomials or the method of undetermined coefficients. We can discard 
the syzygy of which order is more than Y''. It seems to be more efficient to 
determine the coefficients of qi and pi step by step from the low degree. See the 
Example 21 

The membership decision for the order ~<(-w,w,o) needs a tangent cone algo- 
rithm in general, but that for can be performed without the tangent cone 
algorithm. In fact, the membership decision in the step [5] in the Algorithm [2] 
can be performed as follows. 

Algorithm 3 1. Reduce £ G D[s] by {si — Ei \ i = 1, . . . , d}. In other words, 
when £ = J2^aS°', & D, rewrite them to £' := J2 ^aU K' ^ D- 

2. Reduce £' by Ga- Set the resuh to £". 

3. If £" = 0, then £ £ Ha[s] (yes) else £ ^ Ha[s] (no). 

We note that these two steps can be performed in finite steps. If the oper- 
ator £" is zero, then £ G else £ ^ Ha[s]. In fact, when £" ^ 0, then 
in^ ^ (£") ^ F>\s\ ■ {si, . . . , Sd, intu(Gyi)}. It follows from Proposition [5] that 



From Theorems [T] and |4l we obtain the following theorem. 

Theorem 8 The algorithm [3 terminates in finite steps and outputs a correct 
answer. 

Example 3 In the case A ~ (1, 2) and w = (1, 0), the hypergeometric ideal 
Ha[s] = {df — d2,9i + 2^2 — s) has two linearly independent series solutions: 



^Ha[s\. 
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r(^ + i) 



The set of standard monomials determined by Theorem[T]is (1, xidi). We make 
the change of variables yi = — and 2/2 = xi standing for the unimodular matrix 
2 1 \ . . 

^ Q j ' which is the set of the rays of the unimodular cone generated by 

(1,2) and (0,-1). The operators xidi and X2d2 arc transformed to 2yidy^ + 
y2dy2 and —yidy-^ respectively by the transformation. We denote 9y-^ by 9i and 
Oy^ by 62. The hypergeometric ideal and its solutions arc changed as follows: 

R2Ha[s\ = ((^2 + 20i)(02 + 20i-l) + yi0i,02-s), 

Let us derive the Pfaffian system for a direction 9i by following the Algorithm 
[2] and the Remark |4l We put kg — 1 and take a vector F = {l,9i) as a basis. 
We note that (1, 20i + ^2) is a basis from the theorem and then (1, 9i) is a basis. 
Solving the following syzygy equations: 

(911 912) (^J/f^ - (Pii (-s/2) + = 0, 

(911 912) (j^'Si)2jl^ - {Pn P12) ^ ^^^2) + = 0' 

we get the solution 

((711, 912, Pii,Pi2,?'i,(i),i,?'i, (1)^2) = ((0,-4,s(s- l),2(s- 1),0,0), 

(1,0,0,1,0,0), 
(0,0,0,2yi,s,s-l), 
(0,0,yi, 0,1,1)). 

The 3rd and the 4th syzygies are of degree 1 and they should be discarded. 
The first syzygy will be extended to a syzygy of degree 1 in the next step. The 
operator corresponding to the second element is — 0i = G Ha[s]. This 
yields the trivial relation 



mod -ffy!i[s] 



The other operators do not belong to iJ^is]. Next, we put fco = 2 and solve 
above two and the following syzygy equations: 

, X / -.s(s-2)yi/8-s/2\ , .( syi/4+1 \ 2 

'^i^) L(.-2)V/16 + .V4)-(^" ^12) U(,s-2)2;i/8-./2)+"2'(2)^i^i 
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-s(s - 3)?;i/24- (s - l)/2 
^911 9i2j (^(,„i)(,_3)2j^j/48+(s-l)V4 

( \ ( {s- l)yi/12 + l \ 2 

- (m Pi2)(^_(,_i)(,„3)y^/24-(,-l)/2j+"2,(2),2yi-0. 

We get the solution 

(9ll,'7l2,Pll,Pl2, ^1,(1),!, ri_(l)_2, ''2,(2),l' '"2,(2),2) 

((1,0,0,1,0,0,0,0), 

(0, -96, 24s(.s - 1), 24(yi + 4s - 2), 12s, 12(s - 1), 3s(s - 2), (s - l)(s - 3)), 
(0, 48(yi - 4s + 2), -12s(s - l)(yi - 4s + 2), 48(2s - \f , 12s(2s - 1), 12s(2s - 

3s(2s^ - 5s + 4), s - 1, 2s^ - 7s + 9), 
(0, 0, -12y2, 0, 12yi, 3(syi + 4), (s - 1)7/1 + 12)). 

The 4th syzygy is of degree 2 and it should be discarded. The operator corre- 
sponding to the second element is reduced by i?A[s] as follow: 

-966i2 - 24s(s - 1) - 24(yi + 4s - 2)Bx = -24{46l^ + s(s - 1) + (yi + 4s - 2)0^} 

-> -24{46»? + 6'2(6i2 - 1) + {yi + 46I2 - 2)61] 
= -24{(02 + 20i)(02 + 201 - 1) + yi0i} 
^ 0. 

This yields the new relation and we obtain a Pfafhan system for the direction 

y\- 



10\,/1\ /O 1 \/l 



-4/ M6I1 / ~ U(s - 1) yi+4s 



mod 



Example 4 When we use a normal form algorithm by a Grobner basis to obtain 
a Pfafhan system, the matrix Qi is always a diagonal matrix. The degree of 
polynomials in Qi by our Algorithm [5] may be smaller than the diagonal Qi. 

/O 1 1 1\ 

T . •„ . . -.1 IT,.. 100100 

Let us illustrate it with an example, rut A ~ 010010 ^ ^ ~ 

yo 1 1/ 
(si, 1/3, S3, si — 1/2). Note that A ~ 0{P\,2) ~ 1-simplex x 2-simplex. Wc take 

a weight vector w ^ (0,1,1,1,1,0) and y = (ffff , ffff , 2^3, 3^4, 2:5, xe) . The 
basis of series solutions is 

/i = ^"1^1(1/2 -si, -1/3, -S3, 3/2 -si;yi, 2/12/2), 

/2 = 2/''^G2(-si,-S3,l/2,-5/6;-2/i,-2/2), 

/3 - 2/''-^Fi(-S3- 5/6, -SI, -1/3, 1/6; 2/12/2, 2/2), 

where pi are linear polynomials in s^'s, Fi and Gi are the Appell hypergeometric 
function Fi and the hypergeometric function Gi in Horn's list respectively. Here, 
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we apply the Algorithm [5] for the direction 9[ = yidy^. We will denote yi by 
X and 2/2 by y. The degree syzygy returns the pair q = (1, 0, 0),p = (0, 1, 0). 
Here q and p are rows of Qi and P, respectively. The degree 1 syzygy returns 
the pair 

q= (0, + p= (0, --(.T + 3si),Si). 

The degree 2 syzygy returns the pair 

9 = (0, 0, 1 - xy),p - (-^si(2si - 1), i((2si - 2s3 - - 2si + 1), -si) 

where we set the truncation degree N to 5. Note that we get junk syzygies 
conies from lower order syzygies. The space of the degree m syzygies can be 
regarded as a vector space over the base field. Syzygies which comes from lower 
order syzygies can be regarded as a subvector space of the space. These junks 
can be removed by adding an orthogonality condition for the subspace. We 
have used this method to find the second syzygies by removing junk syzygies 
generated by the first one. 

/I \ 

In summary, we have Qi = — a; + 1 x — 1 , of which elements 

\ l~xy j 
have degree 2 at most. On the other hand, if we derive a Pfaffian system with 

/I \ 

a diagonalized Qi, we have Qi = \ (x — l){xy — 1) , of which 

\ l-xy ) 

elements have the degree 3. The Algorithm [5] of utilizing series solutions is 
sometimes slower than the Algorithm[I]in our current implementation. However, 
it often outputs smaller Qi and P,; as we have seen. 

An efficient implementation of the Algorithm [2] for larger problems is not easy. 
This topic together with applications to statistics will be discussed in a forth- 
coming paper. 

Acknowledgments: The authors are grateful to Prof. Maria del Rocio Blanco 
Somolinos for comments on algorithmic resolution of singularities. 
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